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Abstract 

We solve the problem of a Bose or Fermi gas in d-dimensions trapped by S < d mutually perpendicular 
harmonic oscillator potentials. From the grand potential we derive their thermodynamic functions (internal 
energy, specific heat, etc.) as well as a generalized density of states. The Bose gas exhibits Bose-Einstein 
condensation at a nonzero critical temperature T c if and only if d + S > 2, and a jump in the specific heat at 
T c if and only if d + S > 4. Specific heats for both gas types precisely coincide as functions of temperature 
when d + 5 — 2. The trapped system behaves like an ideal free quantum gas in d + S dimensions. For 5 — 
we recover all known thermodynamic properties of ideal quantum gases in d dimensions, while in 3D for S = 
1, 2 and 3 one simulates behavior reminiscent of quantum wells, wires and dots, respectively. 
PACS: 05.30.Fk; 05.30.Jp; 03.75. Fj; 05.70.Ce 

Keywords: Trapped bosons; Bose-Einstein condensation; trapped fermions; quantum wells, wires and dots. 

1 Introduction 

Ultra-cooled bosonic clouds trapped in a harmonic oscillator (HO) external potential mimic the behavior of 
bosons confined by realistic potentials as in opto-magnetic traps in the region of small oscillations. Bose- 
Einstein condensation (BEC) has been now observed with f 7 Rb §, f[Na ||, |Li §, }H {§, ||Rb gjHe § 
and fgK Q neutral bosonic atoms, the upper and lower prefixes being the nuclear mass (number of nucleons in 
the nucleus) and proton numbers, respectively. 
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Trapped quantum gases have been discussed in general by several authors p-p^|. The first calculation of 
the properties of a Bose gas in an isotropic harmonic trap was reported by de Groot et al. ||] ; Bagnato et al. 
^ reported theoretical thermodynamic properties of a Bose gas confined by a generic power-law potential trap; 
Ketterle et al. and Pathria |l3| considered the BEC of a finite number of particles confined in a 3D HO and 
concluded that the thermodynamic-limit approximation is good. For a review of BEC in trapped dilute Bose 
gases see Ref. fl5|| . Trapped Fermi gases have also gained interest as possible precursors of a paired-fermion 



condensate at lower temperatures |16 - Jig], and have been studied experimentally in ultracold fermionic clouds, 
e.g., with |gK neutral atoms in opto-magnetic traps p9[|-p^|. 

On the other hand, the discovery of the quasi-2D superconductors such as cupratc j23)-|2f| or the quasi-lD 
superconductors such as the organo-metallics (or Bechgaard salts) |2^-|2^] have also motivated the study of 
confinement of quantum gases. 

In this paper we describe boson or fermion harmonic-oscillator trapping in order to better understand these 
lower-dimensional structures. Since the system dimensionality modifies the nature of BEC, we seek an exact and 
complete solution in the thermodynamic limit to the non-interacting Bose or Fermi gas problem in (/-dimensions 
constrained by a number 6 of perpendicular HO external potentials. We show that it is possible to map this 
problem into that of a free gas but in a higher dimensionality of d + 5. For example, confinement (29J — J3^] in 3D 
by a ID HO potential collapses the system to a quasi-2D "slab" reminiscent of a quantum well. Confinement 
by a 2D (or 3D) HO potential leads to a wire-like quasi-lD (or dot-like quasi-OD) system. 

In Sec. 2 we calculate the thermodynamic (or grand) potential for the non-interacting Bose or Fermi gas in 
d-dimensions with 5 mutually-perpendicular harmonic traps. In Sec. 3 we deduce the thermodynamic properties 
of these systems and extract a generalized density of states, find their thermodynamic limit, and exhibit their 
mapping to free gases in higher dimensions. In Sec. 4 we specialize to a trapped boson gas and obtain its critical 
BEC temperature, its condensate fraction and specific heat cusp or jump. We also summarize our findings for 
a 3D boson gas trapped by 1, 2 or 3 HO's. In Sec. 5 we specialize to a trapped fermion gas. Sec. 6 contains 
our conclusions. 



2 A cf-dimensional quantum gas trapped by 5 < d HOs 

We consider a (/-dimensional noninteracting boson or fermion gas trapped by 8=1, 2, ...,d mutually-perpendicular 
harmonic-oscillator potentials, the particles otherwise moving freely in the remaining d~ 5 directions. The Hamil- 
tonian for a single boson or fermion of mass m is H — ^2i = iPi /2m + ^muj 2 X)j=d-5+i r j ano - its eigenvalues 
are 

^i^E^^B^ 1 ^) (i) 

where L is the size of the "box" associated with the d — 5 free dimensions and where ni = 0, ±1, ±2, ... while 
Vj = 0, 1, 2, .... We may write the grand potential Q(T, V, fi) in generalized form (see p. 134 of |Q) as 

£l(T, V,(J,) = U — TS — [iN = S a _i^o ~ -S.E. Y ln[l + ae- ,3(e ^^>- Al) ]. (2) 

a ^ — ' 
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Here V = L d ~ s xff is a confinement volume with xo = \Jhjvnuj the oscillator length parameter, U the internal 
energy, T the absolute temperature, S the entropy, /i the chemical potential, N the number of particles, a = — 1 
for bosons, a = 1 for fermions and a — > in the classical case, (5 is the Kronecker delta function, and (3 = 1/ksT . 
In the case of a Bose gas it is convenient to separate out in the sum , the lowest energy state from the excited 
states. We thus defined f^o = — {ksT / a)\o[\ + ae^ 13 ^ 5 / 2 ^^] corresponding to the ground state contribution 
to the grand potential. Using the logarithm expansion ln(l +x) = — YaZi (~ x Y/l valid for |x| < 1, (|j) becomes 

a * — ' * — ' t 

Next, consider only the excited states as the ground state will be treated separately for the boson gas. In the 
continuous limit where H 2 /mL 2 <C kgT and Tiuj -C kgT, the summations over rij and can be approximated 
by integrals, namely E n — * + 1) / d d rii. Thus 



ACT, V, „) = S a - ( 2S + 1 ^ T V f dm e-W^/^nJ 



1 J — OG 



where s is the particle spin, with fermions having s —1/2 and bosons s = 0. The integrals are elementary and 
give 

_ g ~ r_ oe ^( M -*Ba;/2)li 

(M L — mi — • ( 5 ) 



Q_ _i_ I / m r2\( d-d)/2 oo r pfa-snu/l)]! 



The infinite sum is expressible in terms of the polylogarithm function Li a (z) j34j], since 

. > i r 00 , i"- 1 i < az \i 

— aLi a (—az) = — — — / ax — ; = — > — ; . (6) 

v ; r(cr) Jo z- 1 e x + a aj^ l a w 

For a > 1 this reduces to Bose-Einstein (BE) integrals g a {z) when a = — 1 and to Fermi-Dirac (FD) integrals 
f a (z) when a = 1, as defined in Appendix D and E of Ref. |33|, and z = is the fugacity. Using (^|) in (^) 
leaves 

n (T, V, /i) = S^-xClo - ~ ^ (d +g"/ 2+1 ^(<w)/2+i(-a*i) (7) 



where 



_ 2s + 1 , 2 2 ^(d-5)/2 



and 

zi = ze~ f3Shuj/2 = e^Cf- 5 ^/ 2 ). (9) 
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3 Thermodynamic properties 



From (fjj) it is possible to find the thermodynamic properties for a monatomic gas using the relation 

dfl = -SdT - PdV - Ndpi. (10) 

In this representation the grand potential 0(T, V,/i) = — PV is the fundamental relation leading to all the 
thermodynamic properties of the system since 



V - \ — 



1 dT 



dV 



T,fi 



n 
v' 



Using (^) and (|Tl]) the particle number is given by 



N = ~ a( ^dts)/2 L Hd+s)/2{-az 1 ), 



where we used the relation 



Li{d+S)/2{-azi). 



The entropy follows on substituting (^) in the first equation of (|ll|), giving 

S/k B =-[(d+ 5)/2 + 1] a p (d+s) /2 Li {d+s) / 2+ i(-az 1 ) - N In z u 
where we used the number equation ( JT^ ) and the relation 



dLi( d+S )/ 2+1 (-azi) 
dT 



= V (w) L Hd+S)/2(-az 1 ). 
V,n Zl \ U1 / V,n 



Thus (Q) becomes 



C/Ar , [(d + 5)/2 + l]Li {d+s)/2+1 (-az 1 ) 
b/NkB — — t ^ m zi. 



Li( d +s)/2(~azi) 
The internal energy is obtained from (see p. 159 of Ref. |p3[) 

d ( n 



U{T,V) = —kgT 



dT V k R T 



V, z 



Substituting ([?]) here we find that 



and since £1 = — PV then 



U(T,V)=N^- d -±±n, 



PV 



d + 5 



(U-Nhud/2). 



(11) 
(12) 
(13) 

(14) 

(15) 
(16) 

(17) 

(18) 
(19) 



Using (gO and (|12|) the internal energy (|18| ) can be rewritten as 



U(T, V) 
Nk B T 



huiS _^_d + 5 Li( d+ s)/ 2 +i(-az 1 ) 
(—azi) 



2 Li 



(d+S)/2 



The specific heat at constant volume Cy then follows from 



C v = 



^U(T,V) 



N,V 



(20) 



(21) 
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and gives 

CV _ d + 5 
Nk B ~ 2 

where we have used the relation 



,d + S 



Lit 



!+«)/2+i(— d + S Li 



(d+5)/2 



(-azi) 



1 (^1 

zi \ dT 



Li(d+8) /2(- az i) 2 Li^+S)/u-i 

d + S £i(d+5)/2(-a^i) 



(-azi) 



jV.V 



2 £i(d+<5)/2-i(-«^i) 



(22) 



(23) 



which can be extracted from the (vanishing) derivative with respect to T of the number equation (|l2]). Since 

Zi, with g = (d + S)/2 - 1, [d + 6)/2 or 



5 (3(p-6hu/2) > o, (§) then implies that ~aLi a {-az x ) 

(d + <5)/2 + 1. Then @ reduces to 



iVfc B T- 



1 



-21 



1 



(24) 



2(d+<5)/2+l 

which for <5 — 3 gives the classical Dulong-Petit law for T — > oo or zi = 0, while for <5 = we obtain the classical 
limit for ideal gases of bosons or fermions. The first correction to unity in ( pi| ) for d + 6 < 2 is clearly negative 
for a = — 1 (bosons) and positive for a = +1 (fermions), while for d + <5 > 2 it is precisely the opposite. Thus 



we obtain known results obtained for ideal gases for d = [Refs. |35] (bosons), |36 (fermions)]. 

We now recover the results obtained in Refs. p7|-|^9[ dealing with the equivalence of the specific heat 
as a function of T of ideal Bose and Fermi gases in two dimensions. Here this equivalence is obtained more 
generally for d + 5 = 2. If both gases are at the same temperature and have the same number density = uf, 
where = Nb/V is the Bose and uf = Nf/V is the Fermi density, taking d+ S = 2 in ( [l2|) gives 



n B 



A 2 Lix{z\b) 



A 2 Lzi(-zif) 



V 



V =»" (25) 

where as before V was defined just below (j^), z 1B = e t 3 (t J -B-ri^s/2) an j = e ,9(^ F -fiw(5/2) are f U g ac ities 
with fig and /i^? the chemical potentials for bosons and fermions, respectively. Using Landen's relations |3§| ] the 
polylogarithm functions Li a (z) satisfy Li\(x) = —Li\(y) and Li2{x) — —Li2(y) — 1/2 [Lii(y)] 2 , where x — » y 
satisfy the Euler transformation y = —x/(l — x) with a; real < 1. Substituting these relations in (p5|), we obtain 



Zif = (1 - zxb) ■ 



The energy of the Bose gas U(T, V)b taking a = —1 in (20) with d + S = 2, is 



U(T,V) B 
Nk B T 



2 IAx{zib) 



(26) 



(27) 



Substituting (| 



in (^7|) we obtain 

U(T,V) B 
N 



^+^^^ + l/2^Li l{ -z 1F ) 
2 Lii(-zif) 



(28) 

where U(T,V)f is the Fermi gas energy. Substituting (|25| ) in (|28| ) the last term in (|H) is proportional to rip- 
Hence, the energies of the Bose and Fermi gases differ only by a T-independent term and so, from (|2l]), the 
specific heats for boson and fermion gases precisely coincide when d + 5 = 2, or 



[Cy(N > T)] B = [Cy(N,T)} 



F ■ 



(29) 
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3.1 Mapping to higher-d and equivalent mass 



Using (^) and (||), equation (|l2|) can be rewritten as 

N A d+S r dc ^ d+S)/2 ^ _ A d+S r rf£ ( £ -^/2) (d+5)/2 ^ 



T[(d + S)/2]J z^efr + a T ([d + 6} /2) J hujS/2 eK'-n) + a 

deN (e) n (s) , (30) 



where n (e) = [e^ 6- ^ + a] 1 is the BE (a = —1) or FD (a = +1) distribution, and Af (e) is the density of 
states (DOS). Substituting A d+ s from (|J) into ( |30| ) we identify this generalized DOS A/*(e) as 

If 5 = we recover the DOS for a free gas confined in a "box" of sides L 

f mJ 2 \ d/2 pd/2-l 

A6W-P. + i)( 5S ,) rm- ,32 > 

Comparing ( |3l| ) with ( J32] ) in (d + (5)-dimensions 

/mL 2 \ {d+S)/2 p^+^I 2 - 1 

A6(.)-P. + i)( 55P ) rtp+jp)' (33) 



we observe that except for the (negligible) zero-point energy of fko5/2, (|3l| ) and (33) are identical if in (|3l| ) an 
equivalent particle mass m* defined by 

m* = (h/usL 2 ) 25/{d+5) m( d - s )/( d + s ) (34) 

is introduced. Then 

(d+<5)/2-l 

Af(e) = (2s + 1) (m*L 2 /2nh 2 )^/ 2 r{[d + 6] (35) 

In general, therefore, the effect of trapping a quantum gas renormalizes the particle mass m — > m* in accordance 
with ([34]) and increases the dimensionality d — > d + 5 by the number of oscillators. 

3.2 Thermodynamic limit 



Substituting the coefficient A d+ s from (0) into (^) gives 

N = (2s + 1) I ? — ^ - j- / de^ ^J-L 36 

'\2nh 2 J \muj) r([d + *]/2) A w5/2 e^-'O + a 

the volume V being defined just below (j^). The proper thermodynamic limit then holds if N — > oo, L — » oo, 
— > while keeping the ratio N/V — N / L d ~ s x 2 ^ oc Nu> s /L d ~ s = constant. This result was obtained for 
d = 3 and 5 = 3 in Ref. pq| . For a free gas, i.e. 5 = 0, we recover the usual thermodynamic limit, i.e., 
./V — ► co, £ — > oo, while N/L d — constant. 
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4 Trapped bosons 



In this section we study a system of N noninteracting bosons in d dimensions trapped by S (< d) mutually- 
perpendicular harmonic oscillators, and otherwise free in the remaining d — S directions. Let the boson number 
be 

N = No(T) + N k>0 (T) (37) 

were Nq(T) = — (dClo/dfi) T v is the number of bosons in the lowest energy state, with fio defined just below 
(|), while N U>Q (T) is given by (|f) with a = -1. Thus 

N = N (T) + J^ g(d+S)/2 ( Z1 ), (38) 

where from (^) we introduce the Bose function g a {z) which for z = 1 and a > 1 is identical to the Riemann 
Zeta function £(<r). 

Since for T > T c , Nq(T) is negligible compared with N, while for T < T c , Nq(T) is a sizeable fraction of N, 
at T = T c , N (T C ) ~ 0. The critical temperature T c of BEC is found from the condition N k>0 (T c , z\ = 1) ~ N, 
so that ( |38| ) leads to 

N 



2/(d+S) 

(39) 



From (|^) the infinite series g a (l) diverges for a < 1 implying from (|3^) that BEC will occur with critical 
temperature T c ^ if and only if (d + 5)/2 > 1. For 5 — and d = 3 with n = N/L 3 ( |39| ) reduces to the 
familiar formula T c ~ 3.31h 2 n 2 / 3 /mk B of "ordinary" BEC, since o 3/2 (l) = C(3/2) ~ 2.612. On the other hand, 
substituting 6 — 3 and d = 3 in and (0) , we recover the result obtained in |l5[ , 



k B T c ~ 0.94?iw7V 1/3 . (40) 

Enshcr ei aZ., |4(| compared the experimental T c obtained for real traps with the theoretical value ( ^0| ) and 
found good agreement. 

From d38| ) and (|3^) we obtain the condensate fraction, 

N (T)/N = 1 - N k>0 (T)/N (T c ) = 1 - (T /T c ) {d+5)/2 . (41) 

The specific heat follows from ( p2] ) for a = — 1 and from —aLi a (—az) [a, z{\ — g a (z\). We obtain for T > T c 
C v d + S 



^ d + S + 1 ^ ff(d+A')/2+l(^l) _ + £ g(d+(5)/2(^l) 

2 5(d+5)/2(^i) 2 g>(d+5)/2-i(^i) 



iVfc B 2 

while for T < T c , z x = 1 it follows directly from @ and (|l|) that 



(42) 



g v _ ct + '5 ( . <i + (5 | ^(rf+,5)/2 g^d+£)/2+l(l) ^ 



Wfc B 2 2 .g ( d+5)/ 2 (l) 
The specific heat jump at T c is then 

ACV _ C V (T-) - C V (T+) _fd + S\ 2 5(d + 5)/ 2 (l) 



A^s iVfcs V 2 7 g( d +s)/2-i(l)'' 

and, since <7o-(l) diverges for a < 1, will be nonzero if and only if (d + S)/2 > 2. 



(44) 
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The entropy for a = — 1 follows from ([14|). Since —aLi a (^az) [a, z{\ = g a {z\) and using (p2[), in terms of the 
critical temperature T c it becomes 

S/Nk B = [(d + S)/2 + 1] {T /T c f +S)/2 g(d+ * )/2+1 ^ l) - In Zl . (45) 

9(d+5)/2W 

For T < T c , z\ = 1, so that this becomes 

S/Nk B = [(d + S)/2 + 1] [T/T c ) (d+S)/2 g(rf+ ^ )/2 + l( , 1) > 

5(d+5)/ 2 (l) T^Q 

which complies with the third law of thermodynamics. 

For 3D bosons trapped by 1, 2 or 3 harmonic oscillators we summarize our results in Table 1. Since for z% = 1 
the series g a [z\) for a > 1 coincides with C(cr), we require the following values: C(3/2) — 2.612, £(2) ~ 1.645, 
C(5/2) ~ 1.341, C(3) ~ 1.202, C(3/2) ~ 1.127, and ((4) ~ 1.082. In Fig. 1 we show the condensate fraction 
for (5=1,2 and 3. In Fig. 2 are shown their internal energy; specific heat at constant volume (having a jump 
discontinuity if and only if d + S > 4); entropy and chemical potential. 



N /N 
T c 

U/Nk B T 
Cv/Nk B (T < T c ) 
Cv/Nk B (T > T c ) 
AC v /Nk B 
PV 
S/Nk B 



T_\3 

T a 

hull N 11/3 

3 huj i q/ T \3 94(21) 
2 k B T ' 0[ -T c 



12 



C(3) 
C(4) 
CM 

94(^1) _ Q ff3(zi) 



12 



93 (21) 

q03) 
y C(2) 



92(21) 

6.57 



,4/ T 'i3 gj(gi) 
^ItJ C(3) 



lnzi 



i-(£) 5/2 

(2^) 1 / 2 jV xn 
C(3/2) L 



k B T 



2/5 



5^ T ^5/2 97/2(21) 



C(5/2 

35 /'T.V^ ^( 7 / 2 ) 
4 ^T cy l CI372T 

35 97/2(21) _ 25 95/2(21) 

4 95/2(21) 4 93/2(21) 

25 05/2) _ q on 

4 C(3/2) — °- ZU 



Z(U-Nficj) 



Ufj 5/ 



2 97/2(21) 
C(5/2) 



lnzi 



27rjV 2:5 

C(2) T? 



1/2 



7ic£_ , or T \2 93(21) 
:bT ~T ^\T.J C(2) 



1 

2fc B T 1 ~VT, 
2 

6 ' 



/TV C(3)) 

\t c ) ~m 



93(21) 
92(21) 



^ 92(21) 
91(21) 





3(^) 2 W -l^i 



Table 1. Thermodynamic quantities, as 
harmonic oscillators, with the oscillator leni 



defined in text, for a 3D boson gas trapped by 5 = 3, 2,1 
5th parameter xq = {Ti/muj) 1 ^ 2 . 



An ideal Bose gas in d-dimensional space trapped by 8 < d harmonic oscillators has its geometric dimension- 
ality effectively reduced. The BEC temperature expression ( |39| ) for a trapped noninteracting Bose gas shows 
that BEC can occur if and only if (d + S) /2 > 1 as otherwise the term 3(^+5)72(1) diverges forcing T c to vanish. 
Thus BEC is possible in 2D provided S > 1. 

Experiments with dilute boson gases confined in the realistic confining potentials of opto-magnetic traps in 
the region of small oscillations where (BEC) has been observed, can be viewed as a Bose gas in 3D with S = 3. 
Table 2 illustrates some parameters in bosonic vapor systems where BEC has thus far been observed. 
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Boson 
Year/Ref. 


1995 § 


ii Na 
1995 § 


3 Li 

1995 D 


}H 
1998 § 


f 7 Rb 
2000 [§] 


4 He 
2001 § 


41 K 

19 

2001 § 


N 


2 x 10 4 


5 x 10 5 


2 x 10 5 






8 x 10 6 




N 


2 x 10 3 






10 9 


10 4 


5 x 10 5 


10 4 


T c {nK) 


0.17 


2 


0.4 


50 


0.015 


4.7 


0.16 


n (cm -3 ) 


2.5 x 10 12 


1.5 x 10 14 


2 x 10 12 


4.8 x 10 15 


1 x 10 12 


3.8 x 10 13 


6 x 10 11 



Table 2. Some experimental parameters associated with trapped bosonic gases in which BEC 
has been observed to date, TV and No being the number of atoms in the initial cloud and in the 
condensate, respectively; T c the BEC transition temperature; n the boson number density. 



5 Trapped fermions 



Finally, consider a system of N noninteracting fermions in d dimensions trapped by 6 (< d) mutually perpendicu- 



lar harmonic oscillators, and otherwise free in the remaining d—8 directions. Since [e^ e M? 1 )} _|_ j] 



9{E F -e), 



with \i (0) = E F = h 2 k 2 F /2m the Fermi energy, k F being the Fermi wavenumber, we see from ( pO] ) with a = +1 
that 



N 



[2A d+s / (d + S)T {(d + 6)/2\] (E F - hcjS/2) {d+5}/2 



(46) 



[2A d+s / (d + S)T [(d + 8)/2]] E { F d+S)/2 



where in the last step we neglected %u>5/2 compared with E F . The fermion number density with s = 1/2, if 
5 = is obtained ||, ||, substituting (|) in @, as 

n ki 

n = = £ 

L d 2 d - 2 ir d / 2 dT{d/2y 

which reduces to the familiar results n = 2k F /ir, kp/2ir and k F /3w 2 for d = 1, 2 and 3, respectively. 



(47) 



Recalling that —Li a (—z) = f a (z) which are the FD integrals, the internal energy from (GO) can be expressed 



as 



U(T,V) 



P 



fkuS d + 5 f(d+8) /2+i( z i) 



Using (p 



From 



Nk B T 

and the asymptotic expansion for fd/iiz) for T — > (Ref. 
U(T) - NTiuj5/2 (d + S) 



Nk B T F t^o (d + 5 + 2) 

the specific heat as T — > is then 

C V (T) 



2 f(d+8)/2(zi) 

(Ref. H 

(d + 5 



App. B), (fi|) becomes 



1 



2f- ^- 
'\2 \T F 



Nk f 



► (d + S) — 

T^0 6 



T 
TV 



while the entropy S = ^ dT C V (T')/T' becomes 

S/Nk B 



(d + 5f-(l-) 
T^O v 7 6 V^f/ 



(48) 



(49) 



(50) 



(51) 



which again is in agreement with the third law of thermodynamics. These results are reflected in Fig. 3. Table 
3 summarizes results for 3D fermions with 6 = 1, 2, 3 harmonic oscillators. 



!) 



N 

U/Nk B T 
C v /Nk B 
PV/Nk B T 
S/Nk B 



KM" 

3huj 



2k B T 



3 E* 

3/4 (zi) 



1 /4(^1) n /a( z i) 



h(zi) 
4 /-t(zi) _ , 
/ 3 (*i) mzi 



2_g5/2 



!(t^) 1/2 (M- - F 

fig; , 5/ 7 / 2 (zi) 
k B T 2/ 3/2 (zi) 
35/7/2(fl) _ 25/s/2(fl) 
4 /s/2(zi) 4 / 3/2 (zi) 

/7/2(^l) 
/s/2(2l) 
7/7/2(zi) 



2/5/2(2l) 



lnzi 



6 



2fc 

/ 3 (2 



_ 2/3 (Zl) 

T + / 2 (zi) 



/ 2 (zi) ^/l(zi) 

/2(Z1) 



lnzi 



Table 3. Thermodynamic quantities, as defined in text, for a 3D fermion gas trapped by 
5 = 1, 2, 3 harmonic oscillators. 



6 Conclusions 

After constructing the grand potential, thermodynamic properties were determined along with the densities of 
states of ideal boson and fermion gases in d dimensions trapped by S mutually perpendicular harmonic (HO) 
oscillators. Trapping maps the system into a free gas with a new dimensionality increased by the number of 
trapping oscillators, specifically, d — > d + <5, and renormalizes the particle masses m — > m* according to (|34|). 
In particular, we detailed how 3D boson and fermion gases trapped by 1, 2 or 3 mutually-perpendicular HO 
wells map into a free gas in 4, 5 and 6 dimensions, respectively. Also, we found that in a trapped boson gas 
Bose-Einstein condensation with critical temperature T c ^ occurs if and only if d + 5 > 2 so that for S > 1, d 
need not be restricted to d > 2. 
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Figure 1: Condensate fraction for a 3D boson gas trapped by 5 = 1, 2 or 3 harmonic oscillators. 
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T/T c>6 T/T c , 5 

Figure 2: Thermodynamic variables as functions of temperature T, as defined in text, for a 3D boson gas 
trapped by S = 1, 2 or 3 harmonic oscillators. 
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